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The importance of studying the hydrodynamics of the glass-melting tank is substantiated. The main equations 
of a numerical model of the tank hydrodynamics are given. The calculation results of temperature fields and 
streamlines are given for a glass-melting furnace with an output of 300 tons per day. The adequacy of the cal¬ 
culation results as applied to real furnace operating conditions is demonstrated. The developed model of hy¬ 
drodynamics of the melting tank can be used in solving applied problems of designing glass-melting furnaces. 


The evolution of glass-melting furnaces indicates that the 
specific glass melt output is the most critical of numerous pa¬ 
rameters characterizing the engineering efficiency of the fur¬ 
nace. This parameter determines not only the total amount of 
glass produced during the furnace campaign and the thermal 
efficiency of the process, but the furnace cost as well. The 
specific output of furnaces to a large extent depends on the 
melting temperature. The industry considers increasing the 
melting temperature, which is understood as the maximum 
instrumentally monitored temperature of the inner roof sur¬ 
face, as one of the main methods for improving the furnace 
efficiency. At a certain phase of the evolution of glass-melt¬ 
ing furnaces, an increase in melting temperature from 1500 
to 1580°C contributed to raising the glass melt output up to 
2 - 2.5 tons/m 2 per day. However, this is actually the limit¬ 
ing temperature level for the domestic dinas refractory 
(GOST 3910-75) used for the furnace roof. Using high-qua¬ 
lity imported dinas materials makes it possible to increase 
the roof temperature up to 1600°C, and electromelted badde- 
leyite-corandum refractories — to 1640 - 1650°C [1]. How¬ 
ever, in general the possibility of intensifying the furnace 
performance by raising its melting temperature is limited by 
the resistance of refractory materials and is now virtually ex¬ 
hausted. 

Despite the obvious dependence of specific output on 
melting temperature, the latter is a significant, but not the 
only factor determining the furnace efficiency. In our opin¬ 
ion, increasing the specific output of a furnace is an inte¬ 
grated challenge involving a number of technological and 


Ural State Technical University, Ekaterinburg, Russia. 


design problems. A significant factor is the design of the 
melting tank, where the optimization of its geometrical pa¬ 
rameters involves a detailed study of the hydrodynamic regu¬ 
larities of the melt. 

Let us represent the physical model of a glass-melting 
furnace as two zones divided by the glass melt surface and 
separated from the ambient medium by the enclosing refrac¬ 
tory brickwork. In this scheme the technological process 
zone (TPZ, the melting tank) is the main one and the heat 
generation zone (HGZ, the flame space) is the auxiliary zone 
providing for certain conditions in the TPZ. The most impor¬ 
tant prerequisites to obtaining high-quality glass under a pre¬ 
set furnace output include the transfer of a required amount 
of thermal energy into the technological process zone and en¬ 
suring the stay of the melt in the tank for a time required for 
all glass-melting reactions to be completed. The first prere¬ 
quisite is due to the endothermic type of glass-formation re¬ 
actions and the second one is related to the sequence of the 
glass-melting stages and their time dependence. The process 
of heat transfer from the HGZ to the TPZ is the predominant 
one for flame and gas-electric glass-melting furnaces with 
respect to the energy dynamics of the melting tank. The effi¬ 
ciency and the duration of the multistage technological pro¬ 
cess depend both on the amount of thermal energy assimi¬ 
lated by the glass melt and on its distribution in the TPZ vol¬ 
ume. While the first component to a large extent depends on 
the efficiency of external heat transfer, the second one de¬ 
pends on thermophysical properties of the glass batch and 
glass melt and on the hydrodynamics of the melting tank. 
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Fig. 1. Scheme of glass melt migration in the longitudinal section 
of the melting tank of a glass-melting furnace. 


The melting tank hydrodynamic in this case is under¬ 
stood as the set of regularities of heat and mass exchange in 
the tank forming the volumetric temperature field of the 
glass melt. The influence of the external heat transfer param¬ 
eters on the heat and mass transfer within the glass melt vol¬ 
ume are manifested by the temperature profile on the glass 
melt surface. The nature of glass melt flows in the tank de¬ 
pends on several factors. First, it depends on the inhomo¬ 
geneity of the temperature field on the glass melt surface, 
which generates glass melt convection in accordance with 
thermal gravitation convection criteria. Second, the forma¬ 
tion of circulation contours is initiated by the design ele¬ 
ments of the tank impeding natural migration of the melt. 
Finally, the glass melt migration is largely determined by the 
obvious temperature dependence of its viscosity, since the 
tangential stress in the medium forming stream turbulence is 
proportional to the dynamic viscosity coefficient. Thus, we 
recognize a close relationship between the heat and mass 
transfer in the tank and the exterior heat exchange in the 
flame space, as well as the possibility of analyzing them sep¬ 
arately. With such an approach, the results of solving the ex¬ 
terior problem serve as initial data for the interior problem. 
This separation principle has been used in the development 
of a mathematical model for a glass-melting furnace [2, 3]. 

The intensification of heat and mass transfer processes in 
the melting tank achieved by different methods enables one 
to solve other applied problems related to the increasing the 
specific output of the furnace. The point is, primarily, ex¬ 
tending the time of stay of “new” glass melt portions in the 
tank, in particular, by increasing the multiplicity of glass 
melt circulation in the main convection cycles (charging and 
working cycles). This not only decreases the probability of 
substandard glass melt arriving at the working zone, but also 
improves its homogeneity. The increase in the multiplicity of 
the charging convection cycle is especially significant, as 
this has a direct impact on the intensity of the initial melting 
of the batch. The setting of such problems of fluid dynamics 
of the melting tank has become possible due to computerized 
numerical modeling of technological processes. At the same 
time it should be noted that the applied significance of mo¬ 
deling results largely depends on correct formalizing of fluid 
dynamic processes. 

Thermoconvection phenomena in fluids in the general 
case are described by a system of nonlinear partial derivative 
equations. This system includes the equations of the laws of 
conservation of energy, impulse, and mass of material and 


the rheological equation of the state of the medium [4-7]. 
The difficulty of analytically solving fluid dynamics prob¬ 
lems implies the use of various simplifying assumptions. The 
most typical assumption is implied in the problem setting by 
specifying conditions in which the forced motion of the me¬ 
dium is absent and the free motion domain is bounded by 
two parallel walls (top and bottom). Furthermore, the physi¬ 
cal characteristics of the medium are assumed to be constant, 
its density depending only on temperature (and this depend¬ 
ence is taken into account only in the gravity expression), the 
surface temperatures are constant, etc. It is understandable 
that the practical applicability of results obtained with such 
assumptions seems rather conventional. Applied problems of 
designing glass-melting furnaces require a more detailed set¬ 
ting of the problem of the glass melt heat convection neces¬ 
sarily taking into account the exterior (forced) motion of the 
melt through the tank neck. 

Although the temperature field of the glass melt is volu¬ 
metric, at the first stage of the study we can limit ourselves to 
the analysis of the temperature, pressure, and velocity fields 
in the two-dimensional setting (i.e., along the longitudinal 
section of the tank). The validity of such simplification is due 
to the fact that the width of the tank is 3 - 10 larger than its 
depth. Analyzing a stationary (steady) process typical of con¬ 
tinuous furnaces, it is advisable to mathematically describe 
the convection phenomena in the melting tank using the 
“vorticity - stream function” variables. Using the equations 
from [2] the problem for the medium with variable thermo¬ 
physical properties migrating in the area of glass melt motion 
shown in Fig. 1 can be described by the following system of 
equations in dimensionless form: 

vorticity transfer equation: 
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energy equation: 
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where X and 7 are the longitudinal and vertical coordinates 
expressed in fractions of the tank neck height y, (Fig. 1); 
U and V are the longitudinal and vertical components of the 
glass melt velocity expressed in the fractions of its average 
velocity in the tank neck V 0 ; p, p ef , and X ef are the 
dimensionless values of density, effective viscosity, and ther¬ 
mal conductivity of the glass melt; T and c p are the 
dimensionless temperature and specific heat capacity of the 
glass melt; \|/ and ro are the dimensionless values of the 
stream function and flow vorticity; their physical meaning 
ensues from expressions (3) and (4); the dimensionless va¬ 
lues of thermophysical parameters and glass melt tempera¬ 
ture are expressed in fractions of the values of these parame¬ 
ters corresponding to a certain basic temperature T 0 ; Re, Fr, 
and Pe are the Reynolds, Froude, and Peclet similarity num¬ 
bers; the source term R is determined by the change in the 
kinetic energy of the averaged glass melt flow and its effec¬ 
tive viscosity. 

As the area of the glass melt motion has a rather compli¬ 
cated contour, we will use a nonuniform grid (Fig. 2). Let us 
apply the control volume method [8, 9], Having approxi¬ 
mated the differential equations on this grid according to the 
recommendations given in the cited studies, instead of the 
Poisson stream function equation (2) we obtain 

Bp\i/p = Be^e + + ^N^fN + ^S^S + A|/, ptopy (6) 

where the indexes P, E, W, N, and S at the values of the 
stream function \|/ and vorticity ro represent their values in 
the respective nodes of the finite-difference grid (Fig. 2); the 
thermal conductivity values B with indexes P, E, W, N, and S 
include the grid steps along the coordinates X and K as well 
as the dimensionless density of the glass melt at the points e, 
w, n, and s; the coefficient D p is equal to the cross-sec¬ 
tional area of the control volume indicated in Fig. 2 by 
shading. 

The vorticity transfer equation (1) is transformed accord¬ 
ing to the recommendation from [8, 10] using the 
counterflow scheme and S. Patankar’s corrective [10]. As a 
result we obtain a discrete analog of the vorticity equation in 
the following form: 

<n P = Q>, E C0 E + Qo, W (0 W + C(i>,N a N + Q), 5 ra S + An ■ ( 7 ) 

In this case the vorticity conductivity C ro at the nodes P, 
E, W, N, and S in addition to the grid steps include the con¬ 
vective transfer coefficients proportional to the Reynolds 
number and the diffusion transfer coefficients proportional to 
the viscosity p ef . 
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Fig. 2. Pattern of finite-difference grid. 


The glass melt energy equation (5) differs from the vor¬ 
ticity transfer equation only by the presence of a factor at the 
convective terms (specific heat capacity c p ) and the absence 
of the source term. Therefore, the complete finite-difference 
temperature equation can be represented in the following 
form: 

T P = C te T e + C t iyT lr + C tn T n + C ts T s , ( 8 ) 

where the meaning of the thermal energy conductivities C T is 
analogous to C m in Eq. (7). 

To close equation systems (6) - (8), the discrete analog 
of each differential equation in (1) - (5) has to be supple¬ 
mented by algebraic relations approximating the boundary 
conditions at boundary points. The solution of this problem 
requires an explicit representation of the geometry of the 
glass melt motion area, although without yet specifying par¬ 
ticular sizes. 

According to the scheme represented in Fig. 1, the 
boundary conditions on the stream function \\i are strictly de¬ 
fined on the BAEDGHK and IJ contours and in the JK chan¬ 
nel. In the BAEDGHK contour v|/ = 0, whereas the stream 
function distribution in the section JK is determined by the 
glass melt velocity profile. Assuming the glass melt flow di¬ 
rected to the working zone to be laminar, we write 

U(L, Y) = 67(1 - 7), (9) 

where Y = y/y 2 ; L = l/y 2 is the dimensionless melting tank 
length (/ is the tank length). 

Using the first of relations (4) we obtain 

V (L, 7) = \U(L, q)dp = 6^7 2 -^7 3 j = 37 2 - 27 3 , 

where p is the integrating variable. 

Consequently, \|/ = 1.0 on the line IJ. The same value of 
the stream function is preserved on the B ' CF1 line as well. 

According to the mass conservation law, the quantity of 
glass melt arriving at the melting tank on the line BB' should 
be equal to the quantity of glass melt removed at the working 
zone, i.e., the stream function \|/ should vary here from 0 
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Fig. 3. Scheme of boundary parameters. 


to 1.0. The regularity of this variation is determined by the 
kinetics of the physicochemical processes of the formation 
and melting of glass. In the first approximation it can be as¬ 
sumed that the specified processes ensure the required glass 
melt consumption under its parabolic velocity profile, as 
inEq. (9), 
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where X 1 = x l /y 2 (x; is the length of the glass melt area co¬ 
vered by the batch, see Fig. 1); H= h/y 2 is the dimensionless 
height of the glass melt layer (h is the glass melt depth in the 
tank, see Fig. 1). 

The boundary conditions on vorticity co will be specified 
in accordance with D. B. Spalding’s recommendations [8]. 
According to these recommendations, the stream function \|/ 
near the wall changes according to a cubic law, whereas the 
stream vorticity co changes according to a linear law, which 
satisfies the Woods condition [9]: 

_ 3(\\I np ~V|/p ) (% f 


where the indices and the variables correspond to the scheme 
represented in Fig. 3. 

The boundary conditions on vorticity co on the glass melt 
surface can be represented as follows: 

8U dV 

co =-. 

dY dX 


On the line BB' we have U = 0 and, accordingly, 
co = - DV/dX. In accordance with Eq. (10) we can write 


Similarly, V= 0 in the section JK and <o = dU/dY, i.e., 

co (L, 7) = 6(1-27). 

On the line B ' CFI we also have V= 0 and co = dU/dY. 
Considering that this line corresponds to the free surface on 
which the adhesion condition is absent and, consequently, 
there is no vorticity, it is possible in this case to set the zero 
boundary condition. 

Let us consider the temperature boundary conditions. 
The temperature on the glass melt surface is known from the 
solution to the exterior problem, and the heat flux densities 
are known as well. Before solving the problem of the tank 
hydrodynamics, it is necessary just to approximate the zonal 
temperature distribution by a continuous function. This ope¬ 
ration can be performed in any statistical software format. 

On the boundaries AB , AE, IJ and HK the boundary con¬ 
dition of the III kind are specified: 

-Kf^ = k(T-T amb ), 
on 

where n is the dimensionless normal; k is the dimensionless 
heat transfer coefficient; T. is the dimensionless tempera¬ 
ture of the ambient medium. 

Having approximated the derivative in this expression on 
a 2nd order nonuniform grid, we obtain the discrete analogs 
of the boundary conditions on the solid boundaries: 

the left wall (W is the boundary point; dT/dn = - dT/dx) 

Tyy = \.(X E —X W ) Z T E — (X P —X W ) 2 T P + 

k (X F -X p )(X p X w )(X F -X w ) T amb ]/\(X E -X w ) 2 - 
(Xp-X^)- + k (X E —X P )(X P — X W )(X E — X w )\, 

the right wall ( E is the boundary point; dT/dn = dT/dx) 

T e = \{X E -X P f T P - (X E -X w ) 2 T P - 
k (Xp — X W )(X E —X P )(X E — X W )\/[(X E —X P Y — 
(Xp—Xff ,) 2 — k (Xp —X W )(X E —X p )(X E — Xf^)]; 

the melting tank bottom (S is the boundary point, 

dT/dn = - dT/dy) 

T s =[(Y N -Y s ) 2 Tp-(Y P -Y s ) 2 T N + 
k(Y P - Y s )( Y n - Y P )(Y n - Y s )]/[( Y n -7 s ) 2 - 
(Yp— Y s ) 2 + k (Y P — Y S )(Y N — Y p )(Y n — Y s )]. 

One can set adiabatic boundaries as the boundary condi¬ 
tions on the boundaries ED, DG, GH, and JK. In particular, 
when approaching the spillway from the right we obtain 
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Coordinate along the path x , m 

Fig. 4. Glass melt flow in the melting tank (numbers on the curves 
represent normalized stream function). 

The same expression is true for the section JK. When ap¬ 
proaching from the left we have 

T AX E -X w ) 2 T P -jX P -X w ) 2 T E ' 

(X E - X w ) ~(X p - x w ) 

Finally, for the spillway top (line DG ) we have 

T AY n -y s ) 2 t p -(y p -y s ) 2 t n 

s {Y N ~Y s f ~{Y P -Y s f 

Having the complete information on the finite-difference 
setting of the problem, it is possible to select a method for 
solving the discrete analog equations. Due to the substantial 
nonlinearity of the specified equations (all thermophysical 
parameters vary significantly with increasing temperature) 
the only way to solve them is to use the iterative coordinate 
displacement methods. For the same reason the most conve¬ 
nient of the numerous available methods is the consecutive 
displacement method or the Gauss - Zeidel method [11, 12]. 
It is possible to apply various successive over- and under-re¬ 
laxation methods. However, on modem fast personal com¬ 
puters this does not give significant time savings. 

The model of the melting tank fluid dynamics has been 
adapted for a container-glass furnace with an output of 
300 tons per day. The length, width, and depth of the melting 
tank are 13.62, 8.5, and 1.3 m, respectively. The width and 
the height of the tank neck are equal to 0.8 and 0.3 m, and the 
height and the width of the barrier spillway are 0.8 and 
0.4 m, respectively. Other source data are specified in [2], 
The temperature distribution on the glass melt surface was 
taken from the solution to the exterior heat transfer problem. 
The maximum temperature zone is located above the spill¬ 
way installed in the quellpunkt zone. 

To illustrate the calculation results, let us consider the 
distribution of glass melt streamlines (normalized by the 
mass flow rate of the glass melt in the tank neck) and tempe¬ 


y, m 



Fig. 5. Temperature field of glass melt in the melting tank of the 
glass-melting furnace (numbers on the curves denotes glass melt 
temperature, °C). 


ratures in the longitudinal section of the melting tank. The 
data in Figs. 4 and 5 corroborate the quantitative and qualita¬ 
tive adequacy of the obtained results to the generally ac¬ 
cepted concepts of the performance of a glass-melting tank. 
The temperature distribution on the glass melt surface and 
the presence of a spillway contribute to the formation of two 
clearly defined convection contours, the dividing line be¬ 
tween them corresponding to the classical definition of the 
“quellpunkt” notion. 

The intense circulation of the melt is largely due to the 
efficient thermal insulation of the tank brickwork [3] and the 
high specific glass melt output (about 2.6 tons/m 2 per day). 
In comparing the stream function determined by the rate of 
removal of the glass melt to the working zone and a stream 
function equal to 1.0 (in the dimensionless expression), the 
maximum value of the stream function in the tank is 60.66, 
i.e., the circulation multiplicity is more than 60. The intense 
glass melt circulation contributes to homogenizing the tem¬ 
perature in the tank volume. As a consequence, high techno¬ 
logically suitable temperatures are achieved in the bottom 
glass melt layers and in the tank neck. 

Thus, the developed numerical model of hydrodynamics 
of a melting tank combined with the previously described 
model of exterior heat transfer [3] can be used for solving ap¬ 
plied problems in designing or evaluating existing designs of 
glass-melting furnaces. 
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